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ABSTRACT 

We investigate the relationship between black hole mass (Mbh) and Doppler boosted 
emission for BL Lacertae type objects (BL Lacs) detected in the SDSS and FIRST surveys. 
The synthesis of stellar population and bidimensional decomposition methods allows us to 
disentangle the components of the host galaxy from that of the nuclear black hole in their 
optical spectra and images, respectively. We derive estimates of black hole masses via stellar 
velocity dispersion and bulge luminosity. We find that masses delivered by both methods 
are consistent within errors. There is no difference between the black hole mass ranges for 
high-synchrotron peaked BL Lacs (HBL) and low-synchrotron peaked BL Lacs (LBL). A 
correlation between the black-hole mass and radio, optical and X-ray luminosity has been 
found at a high significance level. The optical-continuum emission correlates with the jet 
luminosity as well. Besides, X-ray and radio emission are correlated when HBLs and LBLs 
are considered separately. Results presented in this work: (i) show that the black hole mass 
does not decide the SED shapes of BL Lacs, (ii) confirm that X-ray and optical emission is 
associated to the relativistic jet, and (iii) present evidence of a relation between Mbh an d 
Doppler boosted emission, which among BL Lacs may be understood as a close relation 
between faster jets and more massive black holes. 

Key words: black hole physics - BL Lacertae objects: general - galaxies: active - galaxies: 
jets 



1 INTRODUCTION 

BL Lac type objects and flat-spectrum radio quasars populate the 
class of Active Galactic Nuclei (AGN) denominated as blazars. 
The spectral energy distribution (SED) in BL Lacs usually shows 
a doble-peaked shape. Depending on whether the first peak (syn- 
chrotron contribution) lies in the optical-IR bands or if the peak is 
located at X-ray regime, the BL Lac is classified as low-peaked 
BLLac (LBL) or high-peaked BLLac (HBL), respectively. The 
large majority of AGN detected by the Fermi Large Area Telescope 
in its first 1 1 months of sky survey are blazars and about 40% of 
these sources are classified as BL Lac type objects (Abdo et al. 
2010). As has been shown by recent studies (Lister et al. 2009, 
Savolainen et al. 2010, Tornikoski et al. 2010), 7-ray blazars tend 
to have preferentially higher Doppler factors. Using a sample of 
radio-loud AGN (including mostly blazars), Arshakian et al. (2005) 
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and Valtaoja et al. (2008) have drawn attention to a positive corre- 
lation between black hole mass and Doppler factors. According to 
these two pieces of observational evidence, we may rise the follow- 
ing question: should we expect a population of heavier black holes 
for those blazars detected at high-energies? In order to address this 
question, it is extremely important to obtain reliable estimates of 
the black hole masses in blazars. 

Besides, in order to successfully interpret and model the vari- 
ability in blazars, it is not sufficient to reproduce the observed spec- 
tral energy distributions (SEDs), the variability must occur on a 
physical timescales that is consistent with the chosen model. Thus, 
in order to discriminate between theoretical models and get a true 
estimation of the physical scales in the central engine of blazars, a 
reliable estimation of the black-hole mass must be known. 

When measuring black hole masses in strongly beamed radio- 
loud sources (i.e. blazars) using the gas in the broad line region 
(BLR) and assuming its virial motion (Kaspi 2000, Peterson et al. 
2004), we should be aware of the potential biases in the measure- 
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ments: (i) a possible flat-geometry of the BLR would introduce an 
orientational dependence of the FWHM, (ii) the contamination of 
the optical continuum emission by non-thermal radiation from the 
jets and (iii) non-virialized motions in the BLR (Arshakian et al. 
2010, Leon-Tavares et al. 2010, Ilic et al. 2009). All these potential 
biases make black hole estimations in radio-loud AGN a consider- 
able challenge. 

However, nature is particularly kind and provides us the BL 
Lac type objects, which are those blazars where the host galaxy 
is not outshone by the AGN emission. Thus, enabling us to study 
radio-loud AGN host-galaxy properties and obtain reliable and 
bias-free estimations of the black hole mass via the tight relation 
between the mass of massive black holes and the stellar velocity 
dispersion (Ferrarese et al 2001; Tremaine et al. 2002). 

In this work we perform measurements of the velocity disper- 
sion and bulge magnitudes in a sample of BL Lacs to enable reli- 
able estimations of black hole masses. By comparing our black hole 
estimates to multiwavelegth luminosities and broad-band spectral 
indices, we re-examine the relationship between black hole mass 
and Doppler factors. In the following, we use a ACDM cosmology 
with values within lcr of the WMAP results (Komatsu et al. 2009); 
in particular, H =71 km s _1 Mpc~\ fi m = 0.27, Qa = 0.73. 
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Figure 1. Top panel Spectral synthesis for the BL Lac type object 
J000157.23-103117.3, performed by STARLIGHT (blue line) using 150 
SSP and six power laws to simulate the AGN continuum emission. The ob- 
served spectra (top black line), host galaxy model (green line), AGN com- 
ponent (red line) and residuals (bottom black line) are shown. Bottom panel: 
We show details of the fittings around the absorption features. 



2 THE SAMPLE 

Our BL Lac sample consists of the sub-sample of objects in the 
large sample of BL Lacs from SDSS and FIRST (Plotkin et al. 
2008) for which reliable host galaxy and AGN decomposition has 
been performed in the SDSS spectrum and optical images. The 
Plotkin et al. (2008) sample of BL Lacs was selected jointly from 
SDSS DR5 optical spectroscopy and the FIRST radio imaging. 
Their selection criteria was based on positional matching between 
SDSS and FIRST surveys, and spectral constraints commonly used 
in recent works to classify and characterize BL Lacs, such as no 
exhibition of broad emission lines and a Ca II H/K break measured 
depression C ^ 0.4. 

Within the 256 higher confidence BL Lac candidates in the 
original sample, we have selected 179 objects in the redshift inter- 
val z~ [0.06, 0.50], within which strong stellar absorption lines can 
be measured reliably. As such, it will be in the following identified 
as the original BL Lac sample. In order to separate the spectral con- 
tribution of the host galaxy and AGN, we used the spectral synthe- 
sis of the stellar component method. This allows the star formation 
history, stellar velocity dispersion of the host galaxy and continuum 
spectra of AGN to be modeled and recovered (see section 3.1 for 
details). After a visual inspection of the modeled spectra, 101 out 
of 179 objects have been excluded from the analysis because of the 
poor fitting in the spectral range available, leaving us with 78 BL 
Lacs. We found that the spectra fitting goodness depends mainly 
on the signal-to-noise ratio and the rather moderate contribution of 
AGN optical continuum emission (< 80%). Based on the X-ray to 
radio flux ratio (Padovani and Giommi, 1995), we may classify our 
sample in 45 HBL and 33 LBL. 



3 MEASUREMENTS 

3.1 Stellar velocity dispersion 

The 179 BL Lac spectra from the original sample have been re- 
trieved from SDSS-DR7 and corrected for Galactic extinction us- 
ing the maps of Schlegel, Finkenber & Davis (1998). Then they 



are brought to the rest frame and resampled from 3400 to 9100 A 
in steps of lA with a flux normalization by the median flux in the 
4010-4060 A region. We use the stellar population synthesis code 
STARLIGHlQ to obtain the best fit to an observed spectrum Ox , 
taking into account the corresponding error a bs- The best fit is a 
combination of single stellar populations (SSP) from the evolution- 
ary synthesis models of Bruzual & Chariot (2003) and power-laws 
to represent the AGN continuum emission. 
The code finds the minimum x 2 , 

^ \ (Tabs ) 
X 

where M\ is the model spectrum (SSP and power-laws), obtaining 
the corresponding physical parameters of the modeled spectrum: 
Star formation history, Xj, as a function of a base of SSP models 
normalized at Ao, bj.\, extinction coefficient of predefined extinc- 
tion laws, r^, and velocity dispersion <j* which obeys the relation: 

(N S SP \ 
^2 Xj,b jA r x J <g)G(v*,<T*) (2) 

In order to model the line-of-sight stellar motions, the code 
uses a Gaussian distribution of G centered at the velocity with 
dispersion <r*. We use a base of 150 SSPs plus 6 power laws in the 
form F(A) = 10 20 (A / 4020)' 3 , where /3= -0.5, -1, -1.5, -2, -2.5, -3. 
Each SSP spans six metallicities, Z = 0.005, 0.02, 0.2, 0.4, 1 and 
2.5, Zq, with 25 different ages between 1 Myr and 18 Gyr. Extinc- 
tion in the galaxy is taken into account in the synthesis, assuming 
that it arises from a foreground screen with the extinction law of 
Cardelli et al (1989). 

A detailed description of the STARLIGHT code can be found 
in the publications of the SEAGal collaboration (Cid Fernandes et 
al. 2005; Mateus et al. 2006; Cid-Fernandes et al. 2007; Asari et 
al. 2007). We fit all the wavelength range available in the observed 
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The M BH - Doppler boosted emission relation 3 
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Figure 2. GALFIT decomposition for the BL Lac type object 
J084225. 51+025252. 7. We show (top panel) the 2D original images in the 
g, r and i SDSS-bands, and (bottom panel) residuals after subtracting the 
galaxy model from the original images. 

spectra (3800-9100 A) . To estimate the reliable starlight contribu- 
tion to the optical spectrum we used a uniform weighting to fit all 
the absorption features within the spectral range mentioned above. 
When fitting the observed spectrum, we combined a set of power- 
laws (to account for the AGN continuum luminosity) with stellar 
populations. This allows the reliable simultaneous decomposition 
of the host galaxy and AGN optical continuum emission to be per- 
formed. However, in the case where the host galaxy is under on- 
going star formation we should consider the possibility that AGN 
continuum emission may include a contribution from young stars 
(Cid-Fernandes et al. 2004). 

Table 1 lists the values of the stellar velocity dispersion and 
AGN continuum luminosity (L5100) for those 78 BL Lacs which 
were successfully modeled with STARLIGHT. The uncertainties 
in <r* are the typical errors in the synthesis method, suggested by 
Cid Fernandes (2005) based on the S/N at 4020 A . The spectral 
synthesis result for the BL Lac type object J000157.23- 1031 17.3 is 
shown in Figure 1 , the top-panel shows observed spectra (top black 
line), host galaxy model (green line), AGN component (red line) 
and residuals (bottom black line) are shown. In the bottom-panel , 
we show details of the modeling around the absorption features. 

3.2 Bulge luminosity 

We have modelled the surface brightness profiles of the BL Lac 
host galaxies with the two dimensional image decomposition pro- 
gram GALFIT (Peng et al. 2002). The host galaxy imaging has 
been retrieved from the SDSS photometry images in g, r and i 
bands. Previous studies (Falomo et al. 2000, Scarpa et al. 2000, 
Urry et al. 2000, Nilsson et al. 2003, O'Dowd & Urry 2005 ) have 
shown that the luminosity profile of BL Lac host galaxies is bet- 
ter represented by an r 1//4 law (de Vaucouleurs profile). There- 
fore, we assume a de Vaucouleurs profile to model the host galaxy 
surface brightness in the sample of 78 BL Lac objects selected 
as mentioned in the previous section (3.1). The initial guesses of 
the parameters were obtained using Sextractor (Bertin 1996). Also, 
we used Sextractor to create mask images for the fitting. We se- 
lected bright non-saturated stars in the BL Lacs fields as PSFs. Sky 
background was fitted first and left fixed during model fitting. In 
this way, the total number of parameters is reduced when de Vau- 
couleurs profiles are being computed. 

Every BL Lac galaxy is modeled in the g, r and i SDSS bands, 
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Figure 3. Comparison of black hole mass estimates of 7 1 BL Lac type ob- 
jects from the stellar velocity dispersion and from the R-band bulge magni- 
tude. There is no significant difference in the black hole mass range between 
LBL (open-triangles) and HBL (open-circles) objects. The dashed line re- 
spresents the one-to-one correspondence. 

independently. Then, by comparing the quality of the fittings in 
the 3 bands we decide on the total quality of the fitting. Only the 
galaxies with good fittings in all three bands are considered as part 
of our analysis. In Table 1, we list the 71 BL Lac host galaxies 
best fitted with a simple de Vaucouleurs profile. The bulge magni- 
tude in each band is estimated by using the total magnitude of the 
galaxy and finally the bulge magnitude in R-band (Mr) has been 
computed using the photometric transformation equations found by 
Jester et al. (2005). The final bulge R-band magnitudes have been 
corrected for Galactic extinction (Schlegel et al. 1998) and we ap- 
ply the K-correction from Poggianti (1997). We include in Table 
1 the absolute bulge magnitudes for our best fits. Figure 2 shows 
the results of the structural decomposition for the BL Lac object 
J084225. 5 1+025252.7. 



4 BLACK HOLE MASSES 

Following Tremaine et al. (2002) and McLure and Dunlop (2002), 
the expressions to estimate black hole masses can be expressed in 
the following forms: 

log M sh (ct*) = 4.02(±0.32) log ( —) + 8.13(±0.06) (3) 

\(T0 J 

log Mbh(M r ) = -0.50(±0.02) M R + 2.73(±0.48) (4) 

where ov is the stellar velocity dispersion, o~o = 200 km s _1 and 
Mr is the absolute bulge magnitude in R band (corrected to our 
adopted cosmology). The black hole masses derived from ov and 
Mr are designated as Mbh(&+) and Mbh(Mr), respectively in 
columns 4 and 6 of Table 1 . As can be seen from Figure 3, there is a 
generally good agreement between black hole masses estimated by 
both methods. The average values derived from cr* and Mr are 
(log Mbh) a, = 8.38 ± 0.04 and (log M bh )m r = 8.38 ± 
0.03, respectively, with an average difference of (Alog Mbh) = 
—0.008 ± 0.038 . This result gives us additional confidence in the 
reliability of our black hole mass estimations. Since the relation 
Mbh — suffers of less intrinsic dispersion, is the most accurate 
(Vestergaard 2009) and <r* measurements are not affected by po- 
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28.18 


<25.53 


15.09 


LBL 


J13 1330. 12+020105.9 




254.58zb 13.41 


8.55±0.1 1 


-22.72±0.02 


8.63 ±0 


17 


32.33 


29.06 


< 26.54 


14.20 


LBL 


J 1 32301 .00+04395 1 .3 




280.33 zb 20.30 


8.72±0.15 


-22.28±0.01 


8.41 ±0 


17 


31.67 


28.50 


< 26.90 


18.07 


1 1BL 


J 132617. 70+122958.7 




269.02 zb 16.43 


8.65±0.13 


-21.95±0.01 


8.24 ±0 


6 


31.59 


28.49 


<26.89 


16.32 


11BL 


J133102.91+565541.8 




224.27 zb 22.19 


8.33±0.18 


-21.67±0.03 


8.10±0 


6 


30.88 


28.46 


<26.27 


1 8.66 


1 1BL 


J 133105. 70-002221. 2 




181.24 zbl9.87 


7.96±0.20 


-21.61 ±0.01 


8.07 ±0 


6 


3 1 .21 


28.41 


< 26.00 


16.28 


11BL 


J 134 105. 10+395945.4 




258.51 zb 14.22 


8.58±0.12 


-22. 18zb0.01 


8.36±0 


17 


31.59 


28. 19 


<26.82 


20.06 


1 1BL 


J 134 136.23+55 1437.9 


87GB 133948.0+552941 


222. 53zb 13.45 


8.32±0.12 


-21.90±0.02 


8.22 ±0 


16 


3 1 .('5 


28.65 


<26. 16 


15.39 


LBL 


J 134502.30+5539 1 4.2 




229.72 zb 20.08 


8.37±0.16 


-21.52±0.01 


8.03 ±0 


6 


29.94 


28.39 


24.96 


17.17 


1 1BL 


J 134633.98+244058.4 


NVSS J 1 34634+244 1 00 


209. 73zb 15.24 


8.21 ±0.14 


-22.24±0.01 


8.39±0 


17 


31.37 


28.34 


25 1 6 


12.4 1 


LBL 


J 1353 14.08+3741 13.9 


87GB 135107.4+375518 


289. 96 zb 18.72 


8.78±0.14 


-22.75±0.26 


8.64 ±0 


22 


3 1 .70 


28.76 


<26. 13 


14.79 


LBL 


J 140121. 13+520928.9 




227.32 zb 20.68 


8.35±0.17 


-22.69±0.02 


8.61 ±0 


17 


32.07 


29.06 


<26.87 


16.28 


11BL 


J 140330. 85+36065 1.1 




1 22.47 zb 2 1.26 


7.27±0.32 


-21.20±0.02 


7.87±0 


6 


30.78 


28.50 


25.47 


15.98 


1 1BL 


J 140923.50+593940.7 




258.21 zb 17.87 


8.58±0.14 


-23.45 ±0.01 


9.00±0 


18 


32.27 


29. 1 1 


<27.26 


16.63 


11BL 


J 1 4 1 030. 84+610012. 8 




300.91 zb 23. 53 


8.84±0.16 


-2247 ±0.02 


8.50±0 


17 


31.40 


28.62 


<27.25 


20.25 


1 1BL 


J 142832. 60+42402 1.0 


87GB 142634.5+425353 


260. 35zb 14.53 


8.59±0.12 


-22.14±... 


8.34±0 


17 


3 1 .35 


28.85 


< 27.40 


18.55 


11BL 


J 144248.28+ 120040.2 




310.80zb 13.44 


8.90±0.11 


-22.00±0.01 


8.27±0 


17 


31.77 


29.02 


< 26.98 


16.45 


I1BL 


J 145326.52+545322.4 




21 1.80zb22.71 


8.23±0.20 


-2 145 ±0.01 


8.00 ±0 


6 


29.62 


28. 14 


< 24.78 


17.77 


1 1BL 


J153447.21+371554.5 


87GB 153254.4+372523 


154.11±12.19 


7.67±0.15 


-2 148 ±0.02 


8.01 ±0 


6 


31.15 


28.70 


<25.35 


14.26 


LBL 


J1605 19.04+542059.9 




165.26 ±18.65 


7.80±0.21 


-2 1.02 ±0.03 


7.78±0 


6 


30.90 


28.30 


<26.93 


21.04 


I1BL 


J161541.21+471111.7 




218.22±13.65 


8.28±0.13 


-22.13±... 


8.34±0 


17 


32.09 


28.81 


25.43 


10.63 


LBL 


J 1 6 1 706. 32+4 10647.0 




184.60 zb 9.41 


7.99 ±0.11 


-22.55 ±0.01 


8.54 ±0 


17 


32.16 


29.10 


< 26.71 


14.41 


LBL 


J162839.03+252755.9 




311.76±19.34 


8.91±0.14 


-22.38±0.02 


8.46 ±0 


17 


31.90 


28.72 


< 26.65 


16.28 


11BL 


J163709.50+432600.3 


NVSS J163709+432600 


299.05±11.90 


8.83±0.11 


-23.09±0.01 


8.81±0 


18 


32.29 


29.18 


<25.94 


11.82 


LBL 


J 163726.66+454749.0 


B3 1635+458 


245.11±23.43 


8.49 ±0.1 8 


-22.10±0.01 


8.32±0 


17 


31.28 


28.44 


<26.11 


16.58 


I1BL 


J 1644 19.97+454644.4 


B3 1642+458 


304.30±17.31 


8.86±0.13 


-22.24 ±0.23 


8.39±0.20 


32.27 


28.72 


< 26.63 


14.79 


LBL 


J 165 109. 19+42 1253.4 




175.21 ±18.60 


7.90±0.20 


-2 135 ±0.01 


7.94 ±0 


6 


31.1.0 


28.54 


25.60 


17.82 


LBL 


J171427.40+560 156.0 




268.51±23.41 


8.64±0.17 








31.53 


28.70 


< 26.59 


17.47 


11BL 


J 1729 18.78+525559.2 




175.31±14.80 


7.90±0.16 


-22.85 ±0.36 


8.69 ±0 


25 


31.91 


29.08 


<26.79 


16.58 


I1BL 


J21 1037.75-072412.9 


NVSS J21 1038-072408 


254.84± 17.50 


8.55±0.14 


-2 1.79 ±0.03 


8.17±0 


16 


30.59 


28.47 


25.59 


17.17 


11BL 


J211611. 89-062830.4 




293.34± 13.26 


8.80±0.11 


-22.19±0.01 


8.37±0 


17 


31.61 


28.93 


26.13 


15.39 


LBL 


J221 109.88-002327.4 




273.75±24.05 


8.68±0.17 








32.40 


28.86 


26.70 


14.50 


LBL 


J222944. 18-003426.6 




260.27±12.75 


8.59±0.11 


-22.90 ±0.05 


8.72±0 


18 


30.93 


28.65 


26.82 


20.74 


I1BL 


J235604.01-002353.7 




182.99±20.81 


7.97±0.21 


-2 1.83 ±0.01 


8.19±0 


16 


31.12 


28.68 


< 26.20 


17.47 


11BL 



Table 1. Column description: (1) SDSS name, (2) alias given by radio surveys, (3) stellar velocity dispersions in units (km s _1 ), (4) log of black hole mass 
estimated from stellar velocity dispersion in units of (Mq), (5) R-band bulge absolute magnitude, (6) log of black hole mass estimated from R-band bulge 
luminosity in units of (Mq), (7) log of radio luminosity at 5GHz, (8) log of nuclear optical continuum luminosity at 5100 A , (9) log of X-ray luminosity 
at 1 keV, (10) log of synchrotron peak frequency in (Hz) and (11) the classification of BL Lacs based on the X-ray to radio flux ratio, following Padovani 
& Giommi (1995). Note: All luminosities in this table are in units of erg s Hz -1 and synchrotron peak frequencies in italic font have been taken from 
Nieppola et al. (2006) 



tential beaming effects, in the following we focus our analysis on 
black hole masses estimated via stellar velocity dispersion. 



and LBL objects, confirming results in previous BL Lacs stud- 
ies(Falomo et al. 2002; 2003, Woo et al. 2005) 



Figure 4 shows the distribution of black hole masses esti- 
mated for HBL and LBL objects. As can be seen, the ranges of 
black hole mass in both populations are consistent with each other. 
The Kolmogorov-Smirnov test shows that there is no a significant 
difference between the distributions of black hole mass for HBL 



5 THE BLACK HOLE MASS- LUMINOSITY 
RELATIONSHIP 

In Figure 5 we plot our black hole mass estimations versus mul- 
tiwavelength luminosities. Luminosities at 5 GHz and IkeV have 
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Figure 4. Mgg distribution for 78 BL Lac type objects. There is no a sig- 
nificant difference in the black hole mass distribution between HBL (blue) 
and LBL (red) objects. 



been taken from Plotkin et al. (2008). Disentangling the AGN com- 
ponent from the host galaxy contribution as described in section 3. 1 
allowed us to recover the optical nuclear continuum luminosity at 
5100 A. The solid green lines in the panels of Figure 5 represent 
our fit to the data. Treating Mm as the independent variable, an 
ordinary least-squares fit yields 

r ~, K/r 0.63±0.16 t i\ 

Lsghz oc M BH (5) 

T n,f 0.30±0.11 ,(■-. 

L ' oc M RU (6) 

5100A BH y ' 

T ~, Jl/f 0.48±0.21 tn , 

L lke v OC M BH (7) 

The Spearman rank correlation coefficients listed in Table 2, 
confirm a positive correlation for the above relations at high sig- 
nificance level. These relations are fits for the combined HBL and 
LBL populations. However, the best-fit parameters and correlation 
strengths are slightly changed when HBLs and LBLs are consid- 
ered separately. The fitted relations for LBL and HBL are shown in 
Figure 5 as dotted and dashed lines, respectively. Table 3 lists the 
best-fit parameters for each BL Lac population. It is worth notic- 
ing that the slopes derived for the LBL population are somewhat 
steeper than those derived for the HBL population. 

The correlation between black hole mass and X-ray luminos- 
ity becomes stronger and more significant when computed sepa- 
rately for LBL and HBL populations (see Table 2). This arises nat- 
urally if we bear in mind that HBL and LBL are by definition two 
different populations in the X-ray regime. The Spearman rank cor- 
relation test provides evidence of a positive correlation between 
Mbh and L 5 ghz, L5100 and Li ke v at a confidence level > 98% 
when restricted only to LBL. We note that when considering exclu- 
sively HBLs, the correlations between Mbh and LsGHz,LikeV are 
weaker, though still significant (> 95%). However, the Spearman 
rank correlation test between black hole mass and optical luminos- 
ity for HBL provides only evidence of a trend (at 80% confidence), 
this may betray the presence of a stellar contribution in the opti- 
cal continuum luminosity of HBLs. Whether the stellar component 
differs significantly between LBL and HBL will be addressed in a 
forthcoming paper. As mentioned before, the correlations described 
above have been computed using the black hole mass derived from 
stellar velocity dispersion. However, the best-fit parameters are vir- 
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Figure 5. Correlation between black hole mass and luminosity at 5GHz 
(top panel), 5f00A (middle panel) and IkeV (bottom panel). HBL and 
LBL are shown by open circles and open triangles, respectively. The solid 
line gives the ordinary least-squares fit to the data, whereas fitted relations 
for LBL and HBL separately are shown as dotted and dashed lines, respec- 
tively. See Table 2 for the correlation and best-fit parameters. 



tually unchanged when the estimations of black hole mass via bulge 
magnitude are considered. 

In search of some physical insight into the black hole mass - 
luminosity relation, we explore correlations between X-ray/optical 
and radio emission. For luminosity-luminosity correlations we have 
taken into account the common dependence with the redshift by us- 
ing partial correlation analysis. The results of our partial correlation 
analysis between luminosities are given in Table 2. Radio and opti- 
cal continuum luminosities are shown in the top panel of figure 6, 
where a significant correlation between Lsghz and L5100 is found 
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Table 2. Spearmank's correlation coefficients between black hole mass and 
multiwavelength luminosities for BLLacs in our sample. Partial correlation 
methods have been performed for correlations of the luminosity-luminosity 
type, in order to remove the common dependence with redshift. We consider 
a correlation to be significant if the probability that the correlation is given 
by chance is P ^ 5 X 10 -2 or in other words at a significance level of 95% 
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Table 3. Best fit parameters for the Ms^f-luminosity and luminosity- 






luminosity relations shown in Figure 5 and 6, respectively. 
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isGHz oc M* H 

L 5100A K M lH 

L lkeV oc M| H 

isGHz OC if 10 o 
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0.63 ±0.16 
0.30 ±0.11 
0.48 ±0.22 
1.00 ±0.13 
0.24 ±0.09 



0.70 ±0.21 
0.48 ±0.19 
0.71 ±0.26 
0.86 ±0.14 
0.63 ±0.10 



0.51 ±0.19 
0.18 ± 0.14 
0.48 ± 0.21 
1.02 ± 0.18 
0.60 ±0.10 



at very high confidence level with a correlation coefficient p — 0.6 
and the probability that such correlation is not due to chance is 
> 99.9%) . This suggests that the emission in the two bands has 
a similar origin, more specifically non-thermal radiation from the 
jets. The best fit to the combined sample is shown as a solid line in 
the top panel of Figure 6, whereas dotted and dashed lines represent 
the best fits for the LBL and HBL populations. 

Nevertheless, when attempting to compare the distributions of 
Lsghz and X-ray luminosities, it turns out that we do not find ev- 
idence of significant correlation when the dependence of redshift 
is removed (see bottom panel in Figure 6). However, when HBL 
and LBL are considered separately, the Spearman rank correla- 
tion strength between luminosities becomes higher and significant 
(> 98%). This comes naturally, considering that HBL and LBL 
form two different populations in the X-ray regime. Table 1 lists 
correlation and best-fit parameters for each population. We stress 
that when computing correlations with the X-ray luminosities, we 
use upper limits for those sources non detected in Xrays. It is con- 
ceivable that real values would produce an even higher separation 
between HBL and LBL populations, which would not change the 
significance of correlations computed. 



6 DISCUSSION 

In this work we have measured stellar velocity dispersions and 
bulge luminosities for BL Lac objects identified in the SDSS and 
FIRST surveys, enabling us to estimate masses for the black holes 
in our sample. The range of black hole masses computed from stel- 
lar velocity dispersion is in substantial agreement with those es- 
timated from R-band bulge magnitude, the average difference of 
black hole masses estimated by both methods is {Alog Mbh) = 
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Figure 6. Correlation between radio luminosity at 5GHz and optical contin- 
uum (left panel) and Xray (right panel) emission. The correlation between 
X-rays and radio emission only holds when considering HBL and LBL pop- 
ulations separately. The solid line gives the ordinary least squares fit to the 
data, where fitted relations for LBL and HBL populations are shown as 
dotted and dashed lines, respectively. Best-fit parameters and correlation 
coefficients among luminosities are listed in Table 3. 



-0.008 ± 0.038 . Within the sample of BL Lacs considered, LBL 
and HBL sources are indistinguishable in terms of black hole mass. 

We found that in our sample of BL Lacs, the black hole mass 
correlates at high significance level with radio (5GHz), optical nu- 
clear (5100 A), and X -ray (IkeV) luminosities. Radio luminosity 
at 5GHz correlates tightly with the nuclear optical luminosity and 
with X-ray luminosity when correlation coefficients are computed 
for LBL and HBL, separately, implying that the emission in the 
optical and X-ray bands share a similar origin, more specifically 
non-thermal emission from the jets. 

We note that the strong correlation between black hole mass 
and radio luminosity found in this study is in disagreement with 
the results of Woo et al (2005,) who found no evidence of such 
a correlation using a sample of BL Lacs. The range of black hole 
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Figure 7. The distribution of radio-X-ray spectral index (ot rx ) and syn- 
chrotron peak frequency (v p£a k) for BLLacs which SED has been succes- 
fully modelled in Nieppola et al. 2006. The dashed line gives the ordinary 
least squares fit to the data. 



masses found in this work are similar than those found in Woo et 
al. (2005) and other previous works (Falomo et al. 2003; Dunlop 
et al. 2003). Then, these apparently conflicting results may be at- 
tributable to the literature compilation done by Woo et al. (2005) 
where the measurements have not been obtained in a homogeneous 
way. Whereas in this work all the black hole masses have been es- 
timated in a uniform fashion, data taken with the same instruments. 
Besides, various studies have found a roughly linear correlation be- 
tween black hole mass and radio luminosities (e.g. Falomo et al. 
2003, Dunlop et al. 2003, Arshakian et al. 2005), although con- 
sistent with our study they report a slope higher than ours. This 
apparent disagreement may be related to the mixed population of 
radio sources involved in those studies and also due to the different 
luminosity coverage. The black hole masses used in the previous 
studies are compiled from literature and mainly computed from the 
broad emission lines, which might introduce potential uncertain- 
ties regarding the location, geometry and non-virialization effects 
in the BLR (Arshakian et al. 2010, Leon-Tavares et al. 2010, Ilic 
et al. 2009). The advantage of this work over the previous ones, is 
that our measurements have been performed in a uniform fashion 
and the 78 black hole masses estimated from the stellar velocity 
dispersion are free of the potential uncertainties mentioned above. 

Since it is well established that the synchrotron emission in 
BLLacs (and in blazars in general) is affected by relativistic beam- 
ing effects due to a fast jet aligned close to our line of sight (Bland- 
ford & Konigl 1979), then we must expect the emission in BL Lacs 
to be Doppler boosted. However, we should be aware that luminosi- 
ties in BLLacs might not exclusively depend on Doppler boosting 
effects but on intrinsic physics conditions. 

As a consequence, the correlation Mbh — luminosity found 
in this work, might be understood either as an evidence of observed 
dependence between black hole mass and Doppler boosting effects 
(Arshakian et al. 2005, Valtaoja et al. 2008, Torrealba et al. 2008) or 
as a relationship between black hole mass and intrinsic physics con- 
ditions in each source (i.e. magnetic field strength, jet composition, 
etc). Unfortunately we can not test in situ the correlation between 
black hole mass and Doppler factor with our sample of BL Lacs, 
due to the fact that computation of variability Doppler factors are 
observationally expensive, requiring long-term monitoring. How- 
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Figure 8. Black-hole mass - synchrotron peak frequency plane, where the 
color scale corresponds to the intensity of radio emission, darker colors 
corresponding to brightest radio-jets . The region located at the bottom- 
right of the plane is exclusively populated by the heaviest black-holes 
(M B H > 10 s ' 25 ), brightest radio-jets (L 5 qHz > 10 31 5 ) and the lowest 
synchrotron peak frequencies (u pea k < 10 ie ). Sources contained in this 
beamed-region fullfill the criteria for which broad band emission is domi- 
nated by Doppler boosting effects. 

ever, we may attempt to disentangle the population of BLLacs for 
which luminosity is most likely Doppler boosted, as follows. 

Starting from the hypothesis that emission strength in BL Lacs 
depends primarily on Doppler boosting effects which are charac- 
terized by the Doppler factor (D), then derived expressions (5-7) 
can be seen as L(7J) cx Mbh - Furthermore, it is well established 
that emission from the relativistic beamed jet is responsible for the 
lower energy-peak in the BL Lacs SED and its position on the logv- 
axis is denned as the synchrotron-peak frequency (v pea k). Using a 
sample of 89 blazars, Nieppola et al. (2008) found an inverse de- 
pendence between the Doppler boosting factors and peak frequen- 
cies (D oc v~^ ah )- Such a relation holds for v pea k < 10 16 (Hz), 
above which Nieppola et al. (2008) argued the D~ 1. Based on 
this, we may assume that Doppler boosted emission follow the re- 
lation L(.D) cx v~ k . According to these expressions, we should 
expect the dependence between black-hole mass and synchrotron 
peak frequency to be inversely proportionally (Mbh oc v peak ) for 
BL Lacs where emission is dominated by Doppler boosting effects. 

In order to find whether BLLacs in our sample fulfill the rela- 
tions mentioned above, we compile synchrotron peak frequencies 
estimates for our sample. SEDs for 18 BLLacs contained in this 
study have been constructed and successfully modeled in Nieppola 
et al. (2006). The distribution of peak frequencies with their re- 
spective broad band spectral indexes (chrx) is shown in Figure 7. 
Information about the broad band spectral indexes has been taken 
from Table 6 in Plotkin et al. (2008) . A common least-square fit to 
the data yields, 

a ax = (1.23 ± 0.04) + log(v peak ) x (-0.03 ± 0.01) (8) 

This expression allows us to compute the peak frequencies 
for the rest of our sample. The compiled and computed syn- 
chrotron peak frequencies are listed in Table 1. Figure 8 shows 
the Mbh — v P eak plane where the color scale corresponds to the 
radio luminosity, darker colors corresponding to brightest radio 
jets. We note that there is very well defined region located at the 
bottom-right part of the plane. This region is exclusively populated 
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by the heaviest black-holes(M_g£r > 10 s ' 25 ), brightest radio-jets 
{Lsghz > 10 31 ' 5 ) and the lowest peak frequencies (y pea k < 10 16 ) 
and in the following we refer to it as the beamed region. 32 out of 
78 (40%) BL Lacs in our sample (mostly classified as LBL) are 
enclosed in this beamed region, fullfilling the criteria for which the 
broad-band emission is dominated by Doppler boosting effects. 

Doppler boosting effects in a relativistic jet are often charac- 
terized by the Doppler factor, which is a function of the jet view- 
ing angle and jet speed. Since the jet viewing angle does not play 
a major role in determining the SED shape in BL Lac type objects 
(Padovani & Giommi 1995), we may assume that in our entire sam- 
ple of BLLacs (LBL and HBL), the Doppler boosted luminosities 
will not strongly depend on the orientation angle but mainly on the 
speed of the jet. Hence, the correlation Mbh — luminosity for BL 
Lacs inside the beamed region might be understood as a correspon- 
dence between Mbh and Doppler factor, being the latest mainly 
dependent on the speed of the jet. This can be taken as an evidence 
that more massive black holes produce faster jets in some BL Lac 
type objects. Thus BLLacs outside the beamed region must have 
slower jet speeds. 

Figure 9 shows the Mbh-L$ghz plane where the color-scale 
corresponds to the synchrotron peak frequency. A dashed line has 
been drawn to approximately separate Doppler boosting dominated 
BL Lacs, based on their synchrotron peak frecuency {v pea k ~ 10 16 
Hz). As it can be seen, the separation line coincides with an ap- 
parent dichotomy in radio luminosity, most likely influenced by 
Doppler boosting effects. This implies that all sources above the 
line would shift downwards in Figure 9 if their intrinsic jet lu- 
minosity could be observed. For BL Lacs below the dashed line, 
Doppler boosting should not be a major factor and the observed 
correlation between black hole mass and luminosity is likely to be 
intrinsic. If we consider only sources with v pea k > 10 16 , we are 
still finding that black hole mass correlates (p = 0.3, P = 97%) 
with the radio luminosity, roughly of the form Mbh oc L®q Hz . 
The fact that this slope is shallower than that found for BL Lacs 
with v pea k < 10 16 (Mbh oc L%q Hz ) is merely a consequence 
of Doppler boosting effects on the Mbh — luminosity relation. 
This suggests an additional intrinsic dependence between black 
hole mass and luminosity for sources in which luminosity is not 
Doppler boosting dominated. Whether host-galaxy properties are 
somewhat responsible for the intrinsic luminosity may be gleaned 
from the stellar-population synthesis that comes as a byproduct of 
this study. 



7 SUMMARY 

This work presents black hole mass estimations for a sample of 
78 BL Lac type objects drawn from the SDSS-FIRST sample of 
Plotkin et al. (2008). The main findings of this work are as follows: 

• Within the sample of BL Lacs considered, HBL and LBL are 
indistinguishable in terms of black hole mass. 

• The black hole mass correlates at high significance level with 
radio luminosity at 5GHz, optical continuum luminosity at 5100 
A and X-ray luminosity (see Figure 5). 

• Furthermore, as it can be seen in Figure 6, optical continuum 
and X-ray emission are strongly correlated to the radio luminosity, 
implying that emission in optical and X-ray bands can be associated 
to relativistic jets. 

• When radio luminosity is taken as a third dimension in the 
plane Mbh — v pea k, the sample of BL Lacs separates quite clearly 
in regions where luminosity is dominated by intrinsic processes and 
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log M BH 



Figure 9. Black-hole mass - radio luminosity plane, where the color scale 
corresponds to the synchrotron peak frequency. The dashed line is drawn 
at ^peafe ~ 10 16 Hz, aiming to divide Doppler boosting dominated BL 
Lacs (above the line) from sources where Doppler boosting should not be 
a major factor (below the line) in their luminosity. The correlation between 
Mbh-LsG-STz is present for the overall sample of BLLacs. 

by Doppler beaming effects. The former region is exclusively popu- 
lated by BL Lacs with the heaviest black holes, lowest synchrotron 
peak frequencies and brightest jets. 

• We have found that for a sizeable fraction (40%) of BLLacs 
in our sample, the emission is dominated by relativistic beaming 
effects. For these sources, the correlation between black hole mass 
and luminosity can be understood as an evidence of a correspon- 
dence between black hole mass and Doppler boosting factors. 

• If we assume a typical orientation angle for the population 
of BL Lacs, thus the Doppler boosted emission will not strongly 
depend on the viewing angle but mainly on the speed of the jet. 
Therefore, the correlation between the black hole mass and Doppler 
factor suggested by the results found in this work may imply that 
more massive black holes produce faster jets. 

• For sources which are not Doppler boosted we still find a cor- 
relation between black hole mass and radio luminosity, suggesting 
that it must be intrinsic and more massive black holes produce more 
luminous jets. 

Hopefully, results presented in this work will provide the 
impetus for more observations that look for such correlation among 
blazars. 
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